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We first review how the simple K-matrix unitarized linear SU(2) sigma model can 
Qh, explain the experimental data in the scalar tttt scattering channel of QCD up to 

about 800 MeV. Since it is just a scaled version of the minimal electroweak Higgs 



sector, which is often treated with the same unitarization method, we interpret the 



^ . result as support for this approach in the electroweak model with scaled values of 

tree level Higgs mass up to at least about 2 TeV. We further note that the relevant 
QCD effective Lagrangian which fits the data to still higher energies using the same 
method involves another scalar resonance. This suggests that the method should also 
be applicable to corresponding beyond minimal electroweak models. Nevertheless we 
note that even with one resonance, the minimal K matrix unitarized model behaves 
smoothly for large bare Higgs mass by effectively "integrating out" the Higgs while 
preserving unitarity. With added confidence in this simple approach we make a 
survey of the Higgs sector for the full range of bare Higgs mass. One amusing point 
which emerges is that the characteristic factor of the W-W fusion mechanism for 
Higgs production peaks at the bare mass of the Higgs boson, while the characteristic 
factor for the gluon fusion mechanism peaks near the generally lower physical mass. 
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I. INTRODUCTION 

There has recently been renewed interest in the low energy scalar sector of QCD. 

Especially, many physicists now believe in the existence of the light, broad I = J = 
resonance, sigma in the 500-600 MeV region. The sigma is also very likely the "tip of an 



mixing 



3011 with 



iceberg" which may consist of a nonet of light non qq - type scalars 
another heavier qq - type scalar nonet as well as a glueball. 

A variety of different approaches to meson-meson scattering have been employed to argue 
for a picture of this sort. For our present purpose, first consider an approach to describing 
7T7T, I — J — scattering which is based on a conventional non-linear chiral Lagrangian 
of pseudoscalar and vector fields augmented by scalar fields also introduced in a chiral 
invariant manner Q]. The experimental data from threshold to a bit beyond 1 GeV can 
be fit by computing the tree amplitude from this Lagrangian and making an approximate 
unitarization. The following ingredients are present: i) the "current algebra" contact term, 
ii) the vector meson exchange terms, iii) the unitarized <r(560) exchange terms and iv) the 
unitarized /o(980) exchange terms. Although the p(770) vector meson certainly is a crucial 
feature of low energy physics, it is amusing to note that a fit can be made Q| in which the 
contribution ii) is absent. This results in a somewhat lighter and broader sigma meson, in 
agreement with other approaches which neglect the effect of the p meson. 

For the purpose of checking this strong interaction calculation, the meson-meson scatter- 
ing was also calculated in a general version of the linear SU(3) sigma model, which contains 
both cr(560) and / (980) candidates The procedure adopted was to calculate the tree 
amplitude and then to unitarize, without introducing any new parameters, by identifying 
the tree amplitude as the K-matrix amplitude. This also gave a reasonable fit to the data, 
including the characteristic Ramsauer-Townsend effect which flips the sign of the /o(980) 
resonance contribution. At a deeper and more realistic level of description in the linear 
sigma model framework, one expects two different chiral multiplets - qq as well as qqqq - to 
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mix with each other. A start on this model seems encouraging 

Now, if we restrict attention to the energy range from threshold to about 800 MeV [before 
;he /o(980) becomes important] the data are well fit by the simple SU(2) linear sigma model 
. when unitarized by the K-matrix method If one further restricts attention to the 
energy range from 7T7T threshold to about 450 MeV, the data can be fit by using the non- 
linear SU(2) sigma model 3^|, which contains only pion fields at tree or one loop (chiral 
perturbation theory) level. To get a description of the sigma resonance region by using only 
chiral perturbation theory 3^| would seem to require a prohibitively large number of loops 
(see for example 

The lessons we draw are, first, that the prescription of using the K-matrix unitarized 
SU(2) linear sigma model provides one with a simple explanation of the scalar sector of 
QCD in its non-perturbative low energy region. Secondly, the sigma particle of this model is 
not necessarily just a qq bound state in the underlying fundamental QCD theory. Rather the 
linear SU(2) sigma model seems to be a "robust" framework for describing the spontaneous 
breakdown of SU(2)l x SU(2)r ~ 0(4) to SU{2). There are just two parameters. One 
may be taken to be the vacuum value of the sigma field; this is proportional to the "pion 
decay constant" and sets the low energy scale for the theory. The second may be taken 
to be the "bare" mass of the sigma particle; as its value increases the theory becomes 
more strongly interacting and hence non-perturbative. The K-matrix unitarization gives a 
physically sensible prescription for this non-perturbative regime wherein the model predicts 
the "physical" mass of the sigma to be significantly smaller than the "bare" mass. There is 
in fact a maximum predicted "physical" mass as the "bare" mass is varied. 

Now it is well known that the Higgs sector of the standard electroweak model is identical to 
the SU(2) linear sigma model of mesons. The sigma corresponds to the Higgs boson while the 
7T and 7T° appear eventually, in the Unitary gauge, as the longitudinal components of the W ± 
and Z bosons. There is an important difference of scale in that the vacuum expectation value 
of the Higgs boson field is about 2656 times that of the " QCD sigma" . However, once the bare 
masses of the two theories are scaled to their corresponding VEV's, one can formally treat 
both applications of the same model at once. The practical significance of comparing these 
two applications is enhanced by the Goldstone boson equivalence theorem j^-j^. This 
theorem has the implication that, for energies where My/ j 'E is small, the important physical 
amplitudes involving longitudinal W and Z bosons are given by the corresponding amplitudes 
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of the sub-Lagrangian of the electroweak theory obtained by deleting everything except the 
scalar fields (massless "pions" and massive Higgs boson). Thus, except for rescaling, the 
electroweak amplitudes may be directly given by the pion amplitudes computed for QCD 
with varying bare sigma mass. There is no reason why the scaled Higgs boson mass should 
be the same as the scaled "QCD sigma" mass since, while the two Lagrangians agree, the 
scaled bare mass is a free parameter. 

A simple and perhaps most important point of the present paper is that the same model 
which appears to effectively describe the Higgs sector of the electroweak model can be used, 
with the K matrix prescription to describe the non-perturbative scalar sector of QCD in 
agreement with experiment. In fact, the K-matrix method of unitarization is a popular 
(4]]- (4^ approach to the non- perturbative regime of the electroweak model, although it is 
not easy to rigorously justify. We will see that the QCD sigma meson has about the same 
scaled mass as a Higgs boson of 2 TeV. This provides a practical justification of the use 
of the K-matrix procedure up to that value at least. In fact the successful addition of the 
/o(980) resonance to the low energy sigma using the linear SU(3) sigma model and the same 
unitarization scheme suggests that the range of validity of the electroweak treatment is even 
higher. It also suggests that the same method should also work in models which have more 
than one Higgs meson in the same channel. Furthermore, we will note that, treating the 
unitarized model as a "prescription", even the bare mass going to infinity is a sensible limit. 

A leading experimental question at the moment concerns the actual mass of the Higgs 
particle. Direct experimental search ^| rules out values less than about 115 GeV. Indirect 
observation via a "precision" global analysis of all electroweak data including the effects 
of virtual Higgs boson exchange in loop diagrams actually gives a preferred central value 
of about 80 GeV 4sj|. However there are unexplained "precision" effects like the NuTeV 
experiment on \i neutrino scattering off nucleons j3| and the measurement of the invisible 
Z width at LEP/SLC [45] which raise doubts about what is happening. For example in 
ref J^, the authors suggest the possibility of explaining these two effects by reducing the 
strength of neutrino couplings to the Z boson. In that case, they find a new overall fit which 
prefers rather large Higgs masses, even up to about 1 TeV. This is somewhat speculative 
but does make it timely to revisit the large mass Higgs sector and the consequent need for 
unitarity corrections. 

In the present paper, after giving some notation in section II, we review in section III 
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the fit to low energy I = J = tttt scattering obtained with the K-matrix unitarized SU(2) 
linear sigma model amplitude and its extension to higher energies using an SU(3) linear 
sigma model. This establishes a preferred value for the bare mass of the " QCD sigma" and 
shows it to be in the non-perturbative region of the linear sigma model. The physical mass 
and width, obtained from the pole position in the complex s plane, differ appreciably from 
their "bare" or tree level values. This effect is explored in section IV for a full range of 
sigma "bare" mass values, which can be rescaled to arbitrary choices of bare Higgs mass in 
the electroweak theory. The difference between zero and non-zero pion mass is also shown 
to be qualitatively small. Further, the existence of a maximum value for the physical sigma 
mass is displayed. 

In section V, using the equivalence theorem, we discuss the scattering in the I = J = 
(Higgs) channel of longitudinal gauge bosons for the complete range of bare Higgs masses. 
This provides the characteristic factor in the proposed "W- fusion" mechanism |48[ for Higgs 
production, although here we will mainly regard it as a "thought experiment". Some of this 
material has already been given by other workers ^jj-jl^j]. One feature which is perhaps 
treated in more detail here is the peculiar behavior of the scattering amplitude for large bare 
Higgs mass values. We note that it can be thought of as an evolution to the characteristic 
bare mass = infinity shape. This shape is shown to be simply the K-matrix unitarized 
amplitude of the non linear sigma model. In other words, taking the bare mass to infinity 
in the K-matrix unitarized amplitude effectively "integrates out" the sigma field. We also 
show that the magnitude of the scattering amplitude always peaks at the bare Higgs mass. 
For values greater than about 3 TeV, the width is so great that the peaking is essentially 
unobservable. 



In section VI we discuss the characteristic factor for the "gluon fusion" reaction [491] . which 
is another possible source of Higgs boson production. Unlike the W-fusion reaction, this 
involves final state interactions rather than unit arizat ion. We point out that in our model 
the magnitude of the gluon-fusion factor peaks at the physical Higgs mass rather than at the 
bare Higgs mass found in the W-fusion case. It seems to present an amusing example showing 
how the production mechanism of the Higgs boson can influence its perceived properties. 
Finally some concluding remarks and directions for future work are offered. 
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II. COMPARISON OF NOTATIONS 

First let us make a correspondence between the notations employed for the SU(2)^ x 
SU(2)^ linear sigma model used to model low energy QCD and to model the Higgs sector of 
the minimal standard SU(2)/, x U(l) electroweak theory. In the latter case the Higgs sector 
by itself possesses 0(4) ~ SU(2)i x SU(2)^ symmetry which is explicitly broken when the 
gauge bosons and fermions are taken into account. In the low energy QCD application the 
Lagrangian density is written in terms of the pion and sigma as 

C = ~ (d M 7r • m tt + dpadpa) + a(a 2 + tt 2 ) - b(a 2 + tt 2 )\ (1) 

where the real parameters a and b are both taken positive to insure spontaneous breakdown 
of chiral symmetry. The vacuum value (a) of the a field is related to the pion decay constant 
as 

F w = V2(a). (2) 
(In the low energy QCD model, F„- = 0.131 GeV). The parameter a is given by 

« = \ m lb, (3) 

where m CT b is the tree level value of the sigma mass. Finally the dimensionless parameter b, 
whose value furnishes a criterion for the applicability of perturbation theory, is related to 
other quantities as 

b = ^h- (4) 
8(a) 2 1 ; 

The Higgs sector of the standard model employs the fields 

•=(£) • * t =(*-*°*)- (5) 

These can be rewritten in terms of the tz and a fields by identifying 

* = | " + . | ■ (6) 

V2 

In this language the same Lagrangian, Eq. (0) is written 

C = -0 $t0 $ + - 46($ t $) 2 . (7) 
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Here one employs the notation v = y/2 (4>°) so that, noticing Eq. (J§J) 



(In the electroweak theory, v = 0.246 TeV, about 2656 times the value in the low energy 
QCD case). 



III. UNITARIZED LINEAR SIGMA MODEL 

Now let us consider the SU(2) linear sigma model in Eq. (0) as a "toy model" for the 
scattering of two pions in the s-wave iso-singlet channel. It has been known for many years 
that this gives a good description of the low energy scattering near threshold in the limit 
of large sigma mass. Of course, the non-linear sigma model is more convenient for chiral 
perturbation theory calculations. However we are going to focus on the properties of the 
sigma here and the simple linear sigma model is appropriate for this purpose. 

The I=J=0 partial wave amplitude at tree level is 



1 () 
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(10) 



For generality we have added the effect of a non-zero pion mass which would correspond to 
the addition of a small term linear in a to Eq. The normalization of the amplitude 

Tq(s) is given by its relation to the partial wave S-matrix 



S °( S ) = l + 2zT °( S ). 



(11) 



The formula above is not, of course, new; we are following here the notations of 28(, which 
contains many references. While this tree-level formula works well at threshold it does 
involve large coupling constants and cannot be expected to be a priori reasonable even 
several hundred MeV above threshold. In addition, at the point s = m^ b , the amplitude Eq. 



8 



(JHJ) diverges. A usual solution to this problem is to include a phenomenological width term 
in the denominator by making the replacement: 

— > — = r' 12 

However this standard approach is not a good idea in the present case. As emphasized 
especially by Achasov and Shestakov the replacement Eq. (JT2*j) completely destroys the 
good threshold result. This is readily understandable since the threshold result is well known 
to arise from a nearly complete cancellation between the first and second terms of Eq. Q. 
An advantage of the non-linear sigma model is that the good threshold result is obtained 
directly without need for such a delicate cancellation. However, the pole in the linear sigma 
model can be successfully handled by using, instead of Eq. (|12p. K-matrix regularization 
|3| , which instructs us to adopt the manifestly unitary form 

So{s) ~ i-ipSUM ( } 



Using Eq. (fTT|) we get 



T °( S ) = - [T ° r !l r T r - (14) 



tree 



> S 



Near threshold, where [T$(s)] tree is small, this reduces to Po 1 (s)] iree as desired. Elsewhere 
it provides a unitarization of the theory which is seen to have the general structure of a 
"bubble-sum" . We will adopt this very simple model as a provisional approximation for the 
strong coupling regime of QCD in the I=J=0 channel. 

It seems difficult to rigorously justify this as an effective procedure for the strong coupling 
regime. However we may at least compare with the experimental data on nil scattering j^j] . 
Since F n is known there is only a single parameter - the tree-level mass m^. In Fig. ^ the 
experimental curve for the real part Rq(s) of Tq(s) is plotted up to about y/s = 1.2 GeV. 
Predictions for m ab = 0.5,0.8 and 1.0 GeV are also shown. It is seen that the data up to 
roughly ^/s = 0.8 GeV can be fit when m ab lies in the 0.8 - 1.0 GeV range. Thus, perhaps 
surprisingly, the simple model does a reasonable job of accounting for the low energy, but not 
just the linear threshold region, s-wave tttt scattering data. This circumstance enhances the 
plausibility of the model based on the Lagrangian of Eq. together with the unitarization 
scheme of Eq. (fTIjl . 

There still may be some lingering doubt because the energy region between about 0.8 and 
1.2 GeV is not at all fit by the model. However this is due to the neglect of a second scalar 
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FIG. 1: Comparison with experiment of Real part of the I=J=0 nir scattering amplitude in the 
SU(2) Linear Sigma Model, for m ff j = 0.5 GeV (dots), = 0.8 GeV (dashes) and m a \, = 1 GeV 
(solid). Figure taken from 12811 . Experimental data are extracted from Alekseeva et al (squares) 
and Grayer et al (triangles) 



resonance which is expected to exist in low energy QCD. As shown recently in |28[ if the 
Lagrangian of Eq. (JJJ is "upgraded" to the three-flavor case so that another scalar 

field a' identifiable with the / (980) is contained) the entire region up to about y/s = 1.2 GeV 
can be reasonably fit. This is shown in Fig. [21 in obtaining this fit the value = 0.847 
GeV was selected. Two other parameters, m^f, = 1.300 GeV and a a — a' mixing angle 
were also determined in the fit. Most importantly for our present purpose, exactly the same 
calculational scheme of simply feeding the tree approximation into the unitarization formula 
of Eq. dm was employed. 

In assessing the validity of this approximate model for low energy QCD one should also 
consider the role of the vector mesons. These are known to be important in many low 
energy processes and give the dominant contributions to the "low energy constants" 5f| 
of the chiral perturbation theory expansion. Nevertheless it was found |l4j that, while rho 
meson exchange does make a contribution to low energy s-wave pion pion scattering, its 
inclusion does not qualitatively change the properties of the light a resonance which seems 
crucial to explain the I=J=0 partial wave. More specifically, the effect of the rho meson 
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FIG. 2: Comparison of the best fit for the Real part of the I=J=0 tttt scattering amplitude in the 
non-renormalizable SU(3) Linear Sigma Model with experiment. Figure taken from 28]. Best fit 
value of the bare a meson mass, m a b is 0.85 GeV. 



raises the a mass by about 100 MeV and lowers its width somewhat. 

It also seems worthwhile to remark that this unitarized linear sigma model approach 
gives similar results to a more conventional non-linear sigma model approach [,13] wherein 
the scalar mesons are included explicitly. In order to get the effects of the vector and axial 
vector mesons in the linear sigma model framework one should also add them explicitly in 
a chiral symmetric manner (e.g. j^l). 

IV. SIGMA POLE POSITION 



Let us ask: how non-perturbative is the linear sigma model when it is employed as an 
approximation to low energy QCD? The ordinary criterion for the model to be deep in the 
non-perturbative region is that the dimensionless coupling constant b in Eq. ((!} be very 
much greater than unity. Using Eqs. (J2J) and (j3J) this criterion reads 



\2Fj \2x/2v 



, r , . .. (15) 
\2V2vJ 

Taking m 0.85 GeV to fit experiment, as discussed in the last section, gives a value 
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b = 10.5. Thus it seems fair to say that the theory lies outside the perturbative region 
Nevertheless, the K-matrix unitarization can lead to a result in agreement with experiment. 
In a non-perturbative regime one might expect the physical parameters like the sigma mass 
and width to differ from their "bare" or tree-level values. To study this we look at the 
complex sigma pole position in the partial wave amplitude in Eqs. ()14j) and pOjl : 

t o/„n (ml b - s)a(s) + P(s) 

° [) (ml b -s)[l-ia(s))-i(3( S y 1 ° J 

This is regarded as a function of the complex variable z which agrees with s + ie in the 
physical limit. The pole position zq is then given as the solution of: 

(m 2 6 - zo) [1 - ia(z )) - i(3{z Q ) = 0. (17) 

Note that a(s) remains finite as q 2 = s—Arn^ — > 0, so there are no poles due to the numerator 
ofEq. (Hoi). 

For treating both the low energy QCD as well as the standard electroweak situation it is 
convenient to introduce the scaled quantities : 

m ab m ab 
m = — — 



Z ° = Fl = (18) 

Then, specializing to the unbroken SU{2) x SU{2) case by setting m 2 , = in Eq. (fTU|) 
we may write Eq. (|T7]) explicitly as 



(m 2 - z ){l - i 



— o 

m 
32^ 



-10 + ^ln(l + % 

z m z ' 



}- — = 0. (19) 
J 16tt V ; 



It is easy to solve this equation numerically. In Fig. 3 we show how yRe(,2o) depends 
on the choice of m ab for low energy QCD (solid line). For comparison the situation with 
= 0.137 GeV is also shown (dashed line). We note that the behavior is qualitatively 
similar. 



In Fig. 4, we show how y—Im(z ) depends on the choice of m ab for low energy QCD. 
One sees that the real and imaginary pieces of zq are bounded. Thus we can get an accurate 
analytic approximation to Zq for large m ab by expanding the "log" in Eq. ()19j) . In order to 
get the leading order approximation to Zq it is necessary to keep only three terms: ln(l + x) ~ 
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FIG. 3: Plot of ^/Re(zo) in GeV as a function of m a b in GeV for the QCD application. 

3 . Then, in the SU(2) x SU(2) invariant case with m w = 

352 7T 2 , s 

at large m = . 

The quantity ^jKe(z ) can be interpreted as a physical (renormalized) mass. Now from 
Fig. El and Eq. (|2U|) we note that */Re(zo) decreases monotonically to zero for large rh. 



There is a maximal value of the physical mass, y Re(zo) around m a b = 0.74(0.79) GeV in 
the QCD case with m n = 0(0.137) GeV. This corresponds to the scaled variable rh taking 
the value 5.65 (6.03). Increasing rh beyond this point will decrease the physical mass. Thus 
there are two distinct m's for each value of the physical mass. Low energy QCD selects, in 
the sense of fitting to txti scattering, around 0.85 GeV (with = 0.137 GeV) and the 
corresponding rh around 6.5. The physical mass at that point is around 0.46 GeV. The same 
physical mass would also arise when is about 0.75 GeV. However, as may be seen from 



Fig. HI J — Im(zo) is different for the two cases. At small m a b we may plausibly identify the 
physical width of the sigma as 

For large m„(, we may regard this as a convenient measure of the width, even though 
the simple Breit-Wigner approximation as we will discuss in the next section, is doubtful 
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FIG. 4: Plot of y/—lm(zo) in GeV as a function of in GeV for the QCD application 

there. From Eq. (j2*Uj) we note that, at large m ub [59( the physical sigma mass goes to zero 
proportionally to ^— while the physical measure of the width in Eq. (}2T]) increases as m^. 
Clearly, this limiting behavior is different from a narrow Breit-Wigner resonance. For a 
descriptive understanding of the non-perturbative situation it is probably better to look at 
the unitarized amplitude, Eq. (fT4"j) itself. The real part is illustrated in Fig. [T] and will 
be further discussed in the next section. The mathematical significance of the pole in this 
non-perturbative situation is that the amplitude is approximately given by the sum of a 
complex constant, b a and the pole term: 



Here the residue a a is actually complex. The numerical values of a a and b a are given in Table 
I of Q for the SU(2) linear sigma model case. A check of the accuracy of this approximation 
for the more complicated SU(3) linear sigma model case is illustrated in Fig. 9 of the same 
reference. 




(22) 
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V. APPLICATION TO THE SU(2) x 17(1) ELECTROWEAK MODEL 

Whereas the Lagrangian Eq. (JTJ), characterized by the scale v = GeV, is supposed 
to be an effective theory for calculating the J=0 partial wave amplitudes of low energy 
QCD, the same Lagrangian [written as Eq. (jZJ)] is generally considered to be a portion of 
the minimal electroweak Lagrangian, characterized by the scale v = 0.246 TeV. Of course, 
many people consider this portion to be the least well established aspect of the electroweak 
theory. Furthermore, there is a reasonable probability that the standard model itself is part 
of a still larger theory. Thus it may be appropriate to regard the Higgs sector, Eq. 0, 
unitarized by the K matrix approach, as a kind of effective prescription. An alternative 
"canonical" way to proceed, which has been intensively investigated would be to 

replace it by the non-linear effective chiral Lagrangian, improved by systematically including 
loops and higher derivative terms. However this procedure is expected to be practical only 
to an energy below the Higgs mass (e.g. up to about 0.45 GeV in the analog QCD model, as 
shown in Fig. |2j). Clearly it is desirable to consider a model, like the present one, which has 
the possibility of describing the scattering amplitude around the energy of the Higgs boson 
even if it were to exist in a non-perturbative scenario. 

The discussion of the tttt scattering amplitude Tq in sections 3 and 4 can also be used 
to treat the high energy scattering of the longitudinal components of the W and Z bosons 
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in the electroweak theory by making use of the Goldstone boson equivalence theorem 

This theorem implies that the vector boson scattering amplitudes are related to the 
scattering of the Nambu- Goldstone bosons (with m n = ) of the electroweak theory as 

amp(W£W£ -> WfWZ) = amp(7r + 7r- -> Tr+Tr") + O(^) 

E w 

&mp{W£W L -> Z L Z L ) = amp(7r + 7r- -> vrV) + 0{^-) etc. (23) 

E w 

In the amplitudes on the right hand sides, the value v = = 0.246 TeV should of course 
be used. We will specialize to the Tg(s) partial wave as in Eq. (JTljl and will set = 
in Eq. ([10)1 . This partial wave amplitude will contribute important pieces to the reactions 
in Eq. (J23|) . especially in the vicinity of the Higgs pole. When folded together with the 
appropriate strong interaction pieces the scattering is in principle 4^ measurable from pp 
processes like the one schematically illustrated in Fig. El Of course, competing contributions 
must also be disentangled. 
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FIG. 5: Schematic illustration of pp — > W + W X 

Now, we want to compare the unitarized longitudinal electroweak vector boson scattering 
with the analog ttti scattering. Figs. 3 and 4 show that, in the latter case the effect of m n 
non-zero is qualitatively small. Since < , the effect of non-zero m w is expected 
to be even less significant. Then the Higgs pole positions can be gotten from Figs. 3 and 4 
using the scaled quantities defined in Eq. ()19|1 or the asymptotic formula Eq. (J2U|) . However 
for convenience we display the Higgs pole position in Fig. 6. 

For orientation we recall that the "QCD sigma" has a bare mass m a b of about 850 MeV 
which gives a dimensionless mass m = 6.49 62]. This value of m corresponds to a bare 
Higgs mass value of m a b = 2.26 TeV. At that value, the measure of the physical Higgs mass, 



Ke(zo) would be about 1.1 TeV and y-Im(zo) would be about 1.3 TeV. Evidently the 
QCD sigma is much further in the non-perturbative region than the corresponding range 
of values (< several hundred GeV) which are now usually considered for the Higgs mass. 
Of course there is no reason for the scalar Higgs parameter mo to agree with the QCD 
value. The significance of the present observation is that it increases one's confidence in the 
applicability of the K-matrix unitarization model for Eq. (J7J) to bare Higgs masses in at least 
the 2 TeV region. In other words, the same model with the same scaled parameters agrees 
with experiment there in a different context. Since the present status of the Higgs sector of 
the standard electroweak theory may change if future experiments reveal evidence for new 
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m bare (TeV) m bare (TeV) 

FIG. 6: Square root of modulus of Real and Imaginary parts of pole in unitarized amplitude as a 
function of bare Higgs mass, y/ ReZ po i e m&y be identified as "physical" Higgs mass. Physical width 
may be identified from Eq. Q21|). 

physics just beyond the standard model, it seems worthwhile to have some confidence in an 
approach to a possibly non-perturbative Higgs sector. In particular as the QCD analog two- 
sigma model discussed in section II (see Fig. 2) shows, the K-matrix unitarization method 
can be expected to work in the case of more than one Higgs particle even in non-perturbative 
regions of parameter space. For example, the introduction of the a' resonance which gives 
Fig. 2 would scale to a 2 Higgs model describing physics up to about 3.3 TeV in some 
electroweak theory. 

It may be interesting to give a brief survey of the characteristics of the Higgs particle, 
as "seen" in the s-wave vector boson scattering predicted by the unitarized Higgs model. 
Figure indicates that the physical Higgs mass, \J~Re(z ), approximately agrees with the 
bare Higgs mass m a b until about m ah = 0.5TeV. Thereafter the physical Higgs mass grows 
less quickly |63j until it peaks at about 1.1 TeV with about 2.0 TeV. Afterwards, 
jRe(z ) decreases gradually and vanishes as -^—^ according to Eq. (|2TJj). On the other 
hand, \j— Im(z ) monotonically increases to the saturation value 4t> ^/tt = 1.74 TeV. This 
general pattern of pole position vs. bare Higgs mass has been observed in a variety of 
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At large values of m^, the intuitive meaning of the pole position is not immediately 
clear. It may be more physical to consider the question of the deviation of the predicted 



non-perturbative approaches 
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Higgs resonance shape from that of a pure Breit Wigner resonance. In Fig. 7 are plotted 
|Tq( s )| as a function of ^/s for various values of bare Higgs mass in comparison with the 
pure Breit Wigner shape having the same bare mass and bare width (Tbare = ) • In all 
these plots, the Breit Wigner curve is higher before the peak and lower after the peak. Note 
that the actual bare amplitude in Eqs. Q and (fTUj) includes the effect of crossed channel 
Higgs exchanges and a contact term in addition to the s-channel pole. For a bare Higgs mass 
of 350 GeV there is still not much deviation from the pure Breit Wigner shape. However 
for m a b = 1 TeV the deviation is rather marked. In this case the unitarized model appears 
realistic at lower energies before the peak, which is not surprising since it obeys the low 
energy theorem so is forced to vanish as s — > 0. On the other hand, the simple Breit Wigner 
looks unrealistically high at energies below the peak due to its large width as discussed after 
Eq. ()12|) . Above the peak, however, the unitarized amplitude appears to rise and level off. 
This trend is clarified by the plot for the m CT & = 3 TeV case; there the unitarized model has 
a similar shape to the 1 TeV case below the peak but |Tq(s)| simply saturates to unity for 
higher -y/s. It is also amusing to observe that the peak of |Tq(s)| always occurs at the bare 
Higgs mass, s = vr? ah . This may be seen by noting that, since T° (s) may be expressed in 
terms of the phase shift as exp[z5Q(s)]sin5Q(s), the peak will occur where |5q(s)| = n/2. In 
other words, the peak will occur where Tq(s) is pure imaginary. This is immediately seen 
from Eq (fTT)|) to be the case when s = m 



2 



In order to better understand the large s behavior of the unitarized amplitude it seems 
helpful to examine the real and imaginary parts RefTp (s)] and Ihi[Tq (s)], for various values 
of 771^. Fig. 8 shows these for the same values of the bare Higgs masses as in Fig. 7. Notice 
that the real part RefT^] always vanishes at s = m^ 2 while Iul[Tq] is always unity at that 
point. [These features are evident on inspection of Eq. (JHJ) together with Eq. (JHJ).] While 
this aspect of a simple Breit- Wigner resonance is preserved it is seen that the symmetry 
about the bare mass point gets to be strongly distorted as the bare mass increases. 

One notices that both the real and imaginary parts flatten out at large s for all three 
choices of m a b- This effect may be described analytically by making a large s expansion 
(with = 0) for fixed bare mass m ab of Eqs. Q and (JHJ). [To(s)] tree becomes the 
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(dashed) I=J=0 amplitude for = 0.35, 1, 3 TeV. 
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FIG. 8: Real (solid), Imaginary (dotted) and Absolute value (dashed) of the K-matrix unitarized 
I=J=0 amplitude for m a b = 0.35, 1, 3 TeV. 

By construction, the amplitude is unitary for all s. It is clear that the amplitudes for 
these values of m a ^ vary significantly with energy only up to slightly above m^. The physics 
beyond this point might be filled out if a heavier Higgs meson exists and one would get a 
picture something like Fig. |21 Adding the effects of other interactions involving the Higgs 
meson to [Tq (s)] tree in Eq. (fUj) might also be expected to fill out the flat energy region. 

From the point of view that the K-matrix unitarized model is interpreted as an effective 
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theory (which is unitary for any vn a ^) it is especially interesting to consider the case where 
m„(, gets very large. As a step in this direction, Fig. 9 shows a plot of the amplitude for 
m a \, = 10 TeV. It is seen to be similar to the amplitude for m a b = 3 TeV, except that the 
real part (which does go through zero at ^fs = 10 TeV) saturates to a value which is almost 
indistinguishable from the horizontal axis. In fact this is tending toward a "universal" curve 
- any larger value of will give a very similar shape. This may be seen by expanding 
the amplitude for large m ab while keeping s < m ab . Then [T °(s)] tree « + 0(-£r) so, 

ab 

with K matrix regularization, we get for large m ab 



ReT n u s 



0/ o \ 167TO 2 



0V ' l + (-*-) 2 

1 + I Wttv 2 I 

ImT °(s) 1 —— -3. (25) 

To show the trend towards saturation we have plotted in Fig. EH |^o( s )l f° r the bare 
Higgs masses 1, 3, 10, oo TeV. What is happening should not be surprising. As — > oo 
the tree amplitude is going to nothing but the "current algebra" one, s/(167rt> 2 ). This may 
be obtained directly from the non-linear SU(2) sigma model, which is expected since the 
non-linear model was originally [3^ motivated by taking the sigma bare mass, m ab to be 
very large and eliminating the sigma field by its equation of motion. Thus Eqs. (|23j) just 
represents the K matrix unitarization of the "current algebra" amplitude. Since the result 
is unitary it is not a priori ridiculous to contemplate the possibility that Eqs. (f23J) could 
be a reasonable representation of the physics when m a b describes a far beyond the standard 
model sector of a more fundamental theory. However, all the structure in the scattering 
amplitude is confined to the much lower energy range centered around 4vy/ir ~ 1.74 TeV. 
Of course it may be more likely that the model would represent a unitarization of the theory 
with bare Higgs mass in the 0.1-3 TeV range jfj^]. In any event, the simple m a b 00 
limit nicely explains the evolution of the scattering amplitudes for large bare Higgs mass. 
We should also remark that the non-linear sigma model can be more generally motivated 
directly; there is no need to integrate out the sigma from a linear model. This would lead 
to the alternative "canonical" approach mentioned at the beginning of this section. 
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VI. ADDITIONAL DISCUSSION OF THE ELECTROWEAK MODEL 



We just discussed the Wl — Wl scattering making use of the simple K-matrix unitarization 
and the Equivalence Theorem. This has an application to the "W-fusion" reaction shown 
in Fig. 0D Now we point out that a similar method can be used to take into account strong 



final state interactions in the "gluon-fusion" reaction 
This is an interesting reaction since it is predicted 
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schematically illustrated in Fig. [TTJ 



67 



68, 



69j that gluon fusion will be 



an important source of Higgs production. The t-quark shown running around the triangle 
makes the largest contribution because the quarks, of course, couple to the Higgs boson 
proportionally to their masses. According to the Equivalence Theorem, at high energies this 
Feynman diagram will contain a factor where the 7r's correspond to the Wl's and 

^ ah t 

2 

g a7rn = — ^ appears in the trilinear interaction term ^f^aiz ■ iz obtained from Eqs. (0) and 
(@| . The need for unitarization is signaled, as in the WW scattering case by the fact that 
this diagram has a pole at s — Tn? ab . Now it is necessary to regularize a three point rather 
than a four point amplitude; this is discussed in j^, The WW scattering amplitude 
Pol tree i n t ne P rev i° us section was unitarized by replacing 

. [T °] 



[T °] 



J tree 



tree 



tree 



Portree + ^[To]tree [T(°]tree + 



(26) 




FIG. 11: Gluon Fusion Reaction 
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This has the structure of a bubble sum with the intermediate " Goldstone bosons" on their 
mass shells. For the gluon fusion reaction the final "Goldstone bosons" similarly rescatter 
and we should replace 

m lb ~ s ml b -sl- *[T °] tree 

^-cosage^. (27) 



b 



In the second step we used the relation for the phase shift 5q in the K-matrix unitarization 
scheme: tan^Q = [To] t , where the unitarized I=J=0 S-matrix element is given by S® = 
exp 2i5®. 

It is especially interesting to examine the quantity cos<5q(s) which corresponds to the 
reduction in magnitude of the unitarized gluon fusion amplitude from the one shown in Fig. 
ITTl Using Eq. (jUJ) it is straightforward to find 

cosS° (s) = <b^± _ (2g) 

i(ml b - s) 2 + (a(s) K " s) + Pff 

The numerator cancels the pole at s = m 2 ^ in Eq. (|2T|) so one has the finite result: 

r Qcftt-k r-Ol/ 2 \ 3271V /on\ 

— cos 6%] (s = m ab ) = . 29 

m ah — s 3m a b 

It is clear that the factor cos<5q(s) in Eq. (|27J) effectively replaces the tree level denominator 
{m 2 ab — s) by the magnitude of the quantity in Eq. (|T7Jl . which defines the physical pole 
position in the complex s plane and which we used to identify the physical Higgs mass. 
Note that, as in the W-fusion situation, the regularization is required for any value of 
(not just the large values) in order to give physical meaning to the divergent expression. 
The regularized electroweak factor for the amplitude of Fig. [TTJ 9<yiT \ °° s s ° , is plotted as a 
function of E = y/s in Fig. [12] for the cases m a i> = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 TeV. It is very 
interesting to observe that these graphs clearly show a peaking which is correlated with the 



physical Higgs mass [taken, say, as ijRe(zo)) rather than the bare Higgs mass, m^. At 
rriui, = 0.5 TeV one is still in the region where the physical mass is close to m<j&. Already 
at m„(, = 1.0 TeV the peak has markedly broadened and is located at about 0.89 TeV. At 
rricrb = 2.0 TeV the still broader peak is located near 1.03 TeV, which is about as large as 
the physical mass ever gets. Beyond this, the peak continues to broaden but the location of 
the physical mass goes to smaller s. For example, at m ah = 3.0 TeV, the peak is down to 
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0.88 TeV and is much less pronounced. Going further the curves do not look very different 
from the one with m ab = oo, shown in Fig. ED Here the peaking has disappeared and we 
have the analytic form (recall that g a ^ = ^f 1 )'- 

lim { 9^cos5° (s) 2 + s *-\ (30) 
■m CTb -^oo v m ^ b — s 1<6txv 

This is the analog of Eq. (|25|) and corresponds to the unitarized minimal non-linear sigma 
model. It still takes into account final state interaction effects in direct production of WW 
or ZZ pairs from gg fusion. The non-trivial structure is seen to be confined to the region 
s 1 / 2 less than about 5 TeV. There might be a sense in which such a prescription could be 
appropriate for a situation with an arbitrarily heavy Higgs boson 

It is amusing to compare the effect of the present regularization prescription, Eq. (J27J) 
for the divergence at s = m 2 ub with that of the conventional Breit Wigner prescription used 
for example in 68, 6^, given in Eq. (fT^j) with V = T tree . These are plotted for m ab = 1, 2, 3 
TeV in Fig. El The main observation is that, although the absolute value of the factor 
in Eq. ()12j) flattens out due to the fact that the tree-level Higgs width (to vector bosons) 
increases cubically so the Higgs signal gets lost for larger m a b, the factor in Eq. (j27j) still has 
a peak at lower energies. Also the magnitude of the modified factor suggested in Eq. (|2Tjl is 
larger than the corresponding Breit Wigner factor. Other alternatives to the Breit- Wigner 
prescription for treating the divergence at s = m 2 ab in the gluon fusion Higgs production 
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mechanism in Fig. [n]are momentum- dependent modifications of the Higgs width 
and explicit calculation of radiative corrections 

Of course, the electroweak amplitude factor we have been discussing must be folded 
together with the triangle and gluon part of Fig. ^2 as well as the gluon wavefunctions 
of the initiating particles. Furthermore only the I=J=0 partial wave amplitude has been 
considered. If various partial wave amplitudes are added, the "final state interaction" phase 
factor in Eq. (|27|). exp [«5q(s)] mus t be included too. The value of 5{j(s) may be readily 
obtained from Eq. (J28j) . cos<5q(s) is plotted in Fig. [T51 for representative values of m ab . 
For a detailed practical implementation of this model it would be appropriate to take into 
account the relatively strong coupling of the Higgs boson to the tt channel. This could be 
conveniently accomplished by unitarizing the two channel WlWl — ti scattering matrix. 

On noting that, as we have illustrated (see also Fig. IHJ), there are two bare mass values 
(with different widths) for each physical mass one might wonder if a very light Higgs boson 
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FIG. 12: Unitarized electroweak factor gCT7r T> c ° & , plotted as a function of E = yfs for the cases 
m ab = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 TeV. 



could exist in the strongly interacting mode. Possibly its contribution to "precision elec- 
troweak corrections" would be comparable to those of their light bare mass images. This 
seems interesting, although the corresponding bare masses would be very large (> 10 TeV 
according to Fig. EJ), beyond where the validity of the model has been tested by the QCD 
analog. For bare masses in this region it would seem most reasonable to approximate the 
situation by using the m a b — ► oo case. 
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FIG. 13: Unitarized electroweak factor g<77r7 ^ c ° s ^° , plotted as a function of E = yfs for the cases 

ab 

m ab = 4,oo. 

In this paper we first noted that the K-matrix unitarized linear SU(2) sigma model could 
explain the experimental data in the scalar channel of QCD up to about 800 MeV. Since 
it is just a scaled version of the minimal electroweak Higgs sector, which is often treated 
with the same unitarization method, we concluded that there is support for this approach 
in the electroweak model up to at least Higgs bare mass about 2 TeV. We noted that the 
relevant QCD effective Lagrangian needed to go higher in energy is more complicated than 
the SU(2) linear sigma model and is better approximated by the linear SU(3) sigma model. 
This enabled us to extend the energy range of experimental agreement at the QCD level 
by including another scalar resonance. Similarly in the electroweak theory there are many 
candidates - e.g. larger Higgs sectors, larger gauge groups, supersymmetry, grand unified 
theories, technicolor, string models and recently, symmetry breaking by background chemical 
potentials [z^f - which may give rise to more than one Higgs particle in the same channel. 
We interpreted the better agreement at larger energies in the QCD model as also giving 
support to a similar treatment for a perhaps (to be seen in the future) more realistic Higgs 
sector in the electroweak theory which may be valid at higher energies due to additional 
higher mass resonances. Nevertheless we noted that even with one resonance, the minimal 
K matrix unitarized model behaved smoothly at large bare mass by effectively "integrating 
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FIG. 14: Comparison of modulus of unitarized electroweak factors in Higgs production amplitude: 



(i) 



^7- 

(TO 



for the present regularization scheme (solid) and (ii) prescription of including the tree- 



level Higgs width (dashed) shown in Eq. (|12|). Here the three graphs correspond to m a b = 1,2,3 
TeV. 



out" the Higgs while preserving unitarity. 

With added confidence in this simple approach we made a survey of the Higgs sector for 
the full range of bare Higgs mass. While a lot of work in this area has been done in the 
past we believe that some new points were added. In particular, we have noted that in this 



28 




0.5 1 1.5 2 2.5 3 3.5 4 0.5 1 1.5 2 2.5 3 3.5 4 



ETeV ETeV 




-0.5 - 



-1 - 



0.5 1 1.5 2 2.5 3 3.5 4 

ETeV 



FIG. 15: Plots of cos(5$(s)) for m ab = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0TeV 

scheme the characteristic factor of the W-W fusion mechanism for Higgs production peaks 
at the bare mass of the Higgs boson, while the characteristic factor for the gluon fusion 
mechanism peaks at the generally lower physical mass. 

It should be remarked that while the simplest K-matrix method of unitarization used 
here seems to work reasonably well, it is not at all unique. For example, Eq. ()14|) might be 
replaced by 

° () "l-^([To°] tree + /( S ))' (31) 
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where f(s) is an arbitrary real function. However the fit in Section III is good up to scaled 
: 6. Thus f(s) is expected to be small, at least for this energy range. 



energy, s 



V2- 



A simple modification U in this framework is to take /(.) = Ro[r° Wlli _ Of course 
this is not guaranteed to be more accurate at large s in the non perturbative region. Other 
approaches include using the large N approximation 43^ the Pade approximant method 



I | 0|, the Inverse Amplitude method 



75, 



variational approaches 77] and the 



N/D method |78j. It is beyond the scope of this paper to compare the different approaches. 
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